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Models, where neutrino mass originates from physics at the TeV scale and which are poten- 
tially testable at the LHC, need additional suppression mechanisms to describe the smallness 
of neutrino masses. We consider models in which neutrino mass is generated from the d = 7 
operator LLHuHuH^Hu in the context of SU(5) SUSY-GUTs, where the d = 5 Weinberg op- 
erator can be forbidden by a discrete symmetry. That is, we identify the embeddings in GUT 
multiplets and their consequences for phenomenology and renormalization group evolution. 
We use a specific example to exemplify the challenges. In this case, additional heavy d-quarks 
are predicted, which are constrained by cosmology, in particular, by big bang nucleosynthesis 
and direct searches for heavy nuclei. We show that in the NMSSM extension of the model, 
the discrete symmetry needs to be softly broken, and this soft-breaking can also be the origin 
of deviations from tri-bimaximal mixings. Finally we demonstrate that our example is the 
only tree level decomposition which is consistent with perturbativity up to the GUT scale 
and neutrino mass generation by a leading d = 7 contribution. 



1 Introduction 

The recent measurement of ^13 by Daya Bay [1] and RENO [2] has revealed that the last of the three 
angles of a minimal three- generation model for neutrino physics is relatively large, which opens the 
window for the discovery of leptonic CP violation. This result exerts pressure on flavor symmetry models 
(such as based on the permutation groups and 54), as they are generically constructed in a way 
to give at leading order tri-bimaximal (or bi-maximal) mixing - thus implying ^13 = (see [3] for an 
extended review). Corrections from the charged sectors or next-to-leading contributions to the neutrino 
mass matrix are generally invoked to correct such patterns and make the models compatible with the 
experimental data. 

With the recent discovery of a bosonic resonance [4, 5], showing all the characteristics of the SM Higgs 
boson, a long search might soon come to a successful end. In contrast no conclusive hints for physics 
beyond the Standard Model (SM) have been found [6-15]. Therefore, large areas of the parameter space 
of the simplest SUSY models are excluded now and the allowed mass spectra as well as the best fit values 
to the data are pushed to higher and higher values, see, e.g. [16, 17] and Refs. therein. This has lead to an 
increasing interest in the study of SUSY models which provide new features. For instance, models with 
light stops and higgsinos can more easily explain the recent findings within the minimal-supersymmetric 
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standard model (MSSM), see e.g. [18-21], while models with broken i?-parity might be able to hide much 
better at the LHC [22, 23]. Moreover, the MSSM extension with a singlet Higgs boson [24-29] or an 
additional gauge group can easier explain a Higgs boson with a mass of 126 GeV than the MSSM [30-36]. 
Most of these models will be tested at the LHC after its energy upgrade to 14 TeV. 

Generic models to describe the smallness of neutrino mass, such as the standard type-I seesaw mecha- 
nism [37-40], imply new physics close to the GUT scale. In these seesaw models, one can integrate out 
the heavy mediators in order to obtain the famous d = 5 Weinberg operator [41], 

Cif^-;^(L-iT'H)imr'L) (1) 
Anp 

which leads, after Electroweak Symmetry Breaking (EWSB), to Majorana masses m,^ ~ w^/Anp for the 
neutrinos. Here L and H stand for the SM lepton and Higgs doublets, respectively, v for the Higgs vev, 
and Anp for the new physics scale. For order one couplings, Anp points to the GUT scale, which means 
that these models are not accessible by experiments. In addition, they contain more parameters than 
there are observables in the neutrino sector. 

It is therefore potentially attractive to test neutrino mass models with new particles at the TeV scale, 
since current and future LHC upgrades can directly discover or constrain these models. In such models, 
the smallness of the neutrino mass requires (at least) an additional suppression mechanism: Examples 
are radiative mass generation, where loop suppression factors enter, see Refs. [42-47] for systematic 
approaches, or models with a small lepton number violating contribution, such as the inverse seesaw 
mechanism [48] or SUSY with R-parity violation where both aspects enter, see for example [49-58]. 

In this paper, we focus on another possibility, which has been recently drawing some attention: If the 
operator in Eq. (1) is forbidden, such as by a discrete symmetry, a higher dimensional operator may 
dominate [43, 59-70]. Such models can also be combined with supersymmetry, see for example [61, 71]. 
The simplest possibilities of operators in the MSSM, leading to neutrino mass, might be 



1 

A 



r.t^^''^' = -iz^iLLMMiH^HuT , n = 1,2,3,..., (2) 



whereas in the NMSSM, 



NP 



Ctr^' = -^{LLH^Hu){Sr , n= 1,2,3,... (3) 



Anp 



or mixed combinations with Eq. (2) are possible, see Ref. [71] for a detailed discussion. In all these cases, 
neutrino mass is suppressed by a higher power of Anp, which means that Anp may be potentially the 
TeV scale. 

The simplest possibility which works both in the MSSM and the NMSSM is the d = 7 operator 

^if = -A- {LLHMiHdHu) . (4) 

^^NP 

This operator is interesting for two reasons: First of all, it already implies a substantial suppression 
of neutrino mass, more than the d = 6 operator in Eq. (3) (for n = 1). And second, in the NMSSM 
framework the concept of neutrino mass generation can be extended to higher [d > 7) dimensions, c/., 
Eq. (2), whereas Eq. (3) as leading contribution to neutrino mass is in fact limited to n < 2, and operators 
with mixed combinations of S and {HdHu) cannot exceed d = 7 as well.^ Therefore, we focus on Eq. (4) 
in this study. 

In Ref. [71], we have studied a particular example for the generation of this operator which leads poten- 
tially to interesting signals at the LHC. In this paper we go a step further and investigate how potential 
GUT completions of such models look like and in which SU (5) representation the corresponding media- 
tors would reside in. This is motivated by the fact that some of the newly postulated particles are charged 



^The reason for that in the NMSSM is that the SHuHd and S"^ terms have to be allowed by the discrete symmetry, 
which implies that certain lower dimensional operators leading to neutrino mass cannot be forbidden [71]. 
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Figure 1: Example for neutrino mass chosen in section 2. Here N and N' are fermion singlets, and ^ and C 
are fermion doublets. 



under the SM gauge group and, thus, would destroy the successful MSSM prediction for gauge coupling 
unification [72-78]. This is similar to the case of the seesaw type II and type III mechanisms in which 
the corresponding SU (2) triplets have to be included in complete SU (5) multiplets leading to interesting 
effects for collider physics and the abundance of the relic dark matter density [79-88]. The requirement 
to obtain gauge coupling unification implies that not only the mediators for the higher dimensional neu- 
trino mass operators have masses at the TeV scale, but also the members of the corresponding SU{5) 
multiplets. However, renormalization group effects can lead to a sizable mass splitting in particular for 
those particles charged under SU{3)c- 

The plan of this paper is as follows: we start with the model presented in Ref. [71] and discuss its GUT 
completion considering both, the MSSM Higgs sector as well as the NMSSM Higgs sector. As we will 
see, we have to embed heavy SU{2) leptons, required by the model, into SU{5) five-plet representations, 
which implies the prediction of new heavy d-quarks at the TeV-scale. In order to avoid the /i-problem of 
the MSSM, we discuss the NMSSM implementation in section 3 and its challenges, as well as we propose 
an origin of deviations from tri-bimaximal mixings in this approach. In section 4, we review the LHC 
and cosmological constraints on the heavy d-quarks, and propose fast enough decay modes compatible 
with the model. In section 5 we systematically discuss the GUT completions of the possibilities for 
neutrino mass from the higher dimensional operator suggested in [71]. Here we show that several of them 
would imply non-perturbative gauge couplings already significantly below the GUT scale. We constrain 
ourselves to the SU{5) case as it is also a subgroup of larger groups such as 50(10) or Eq. Finally in 
section 6 we summarize and draw our conclusions. 



2 A low energy model for a realization of a d = 7 neutrino mass operator 
and its GUT completion 

In Ref. [42], it has been demonstrated that at tree-level, there are exactly three decompositions of the 
Weinberg operator, which are known as type I, type II, and type III seesaw mechanisms, respectively. 
Similarly, the d = 7 operator in Eq. (4) can be decomposed systematically at tree-level, as discussed 
in Ref. [65] for the Two Higgs Doublet Model and in Ref. [71] for SUSY. As a side remark, note that 
supersymmetry constrains the possible decompositions of these operators in contrast to the Weinberg 
operator, where all three possibilities have also a supersymmetric realization. 

Here we will first recall a particular decomposition of a d = 7 operator, which has been already 
presented in Ref. [71], and discuss its GUT completion to exemplify the challenges one encounters. The 
tree- level diagram for the leading contribution to neutrino mass is shown in Fig. 1. The corresponding 
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superpotential is, at the electroweak scale, given by 



W = W^^ssM + YnNL ■ Hu - f^iN'i ■ Hd + K2N'i' ■ + tumNN' + m^i ■ |' , (5) 

where the mediators N and N' are SM singlets and ^ and ^' are SU{2)l doublets carrying hypercharge 
1/2 and —1/2, respectively. After electroweak symmetry breaking the mass matrix for the neutral weakly 
interacting fermions reads schematically in the basis = {lyjNjN',^^,^'^) as 



MO 






YnVu 











YnVu 





IJlN 











rriN 




























K2Vu 








(6) 



\ K2Vu ni^ / 



where we have neglected all flavor indices. By integrating out the heavy fields one obtains an effective 
mass matrix for the three SM neutrinos 

= vlvdY^ , (7) 

where the couplings carry a flavor index. As has been demonstrated in Ref. [71], one successfully explains 
neutrino masses and mixings within this model by an appropriate choice of fermion generations and 
coupling matrices. In addition, one can easily see from Eq. (6) that the neutral fermion mass matrix 
reduces to an inverse seesaw form if the and ^''^ are integrated out, where the /i-term (the 3-3-element 
in the inverse seesaw mass matrix) is suppressed by m^. 

Choosing the masses of the mediator fields at the TeV scale implies that the additional couplings are 
0(10"^), which is comparable to the SM Yukawa couplings. Accordingly, this model can be tested at 
the LHC. The SM singlet fields and N' are only produced in small amounts due to the smallness of 
the Yukawa couplings. The SU(2) doublets ^ and however, can be produced in Drell-Yan processes, 
similarly to charginos and neutralinos, with a cross-section of up to 0(10^ fb). These particles will then 
decay into vector bosons and leptons which is induced by the mixing between the neutral fermions as 
described by the mass matrix in Eq. (6). The smallness of the mixing between the heavy and the light 
neutrinos implies that one has to expect sizable decay length of up to several millimeters. Under favorable 
conditions one can even establish connections between neutrino physics and lepton number violating LHC 
signals arising from the production of two same sign charged leptons and several VF-bosons. 

The presence of the ^ and ^' supermultiplets modifies the running of the gauge couplings such that 
unification would not occur anymore. Therefore, one has to embed them into a complete SU{5) repre- 
sentation and, as they carry the same quantum numbers as the MSSM Higgs bosons the natural choices 
are 5^ and 5^'. The minimal field content is given by 

-KT). ^.-(n^ ^-(t')' --(t)' --(t)' 

and the fermionic singlets N and N' . This implies that at the electroweak scale one has not only the 
new leptonic fields and ^, but also additional heavy quarks d'^ and d" stemming from 5^/ and 5^, 
respectively. We will discuss their phenomenology in section 4. For completeness we recall that the fields 
A'" and A^' are (fermionic) gauge singlets, similar to the NMSSM singlet scalar field S. 

The most general SU{5) invariant superpotential containing the MSSM fields, 5^', 5^, A^ and A^' reads 

as 

W = yiN5^H^ + y2N5^, H5 + y3N5MH5 + 

y[ , N' 5^ % + y'2 N' 5^, H5 + y'^ N' 5m ^5 + 

m^/ 5m 5^ + 5^/ 5^ + itinN'N + uimnNN + uim'N'N'N' + 

yd 5m 10 % + y'd 5^' 10 + y„ 10 10 - fiH-^H^ . (9) 
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Table 1: Possible Z3 assignments to forbid the Weinberg operator. 



Multiplet 


5a/ H5 i?5 N N' 5^ 5^/ 10 


Z3 charge 


1 1 1 1 2 1 



If both ys and y'^ were present, then one would generate the Weinberg operator which in general would 
be the dominant source for neutrino masses. However, if there were a discrete flavor symmetry operating 
at the GUT scale, then this operator could be forbidden, implying that the d = 7 operator in Eq. (4) is 
the dominant contribution to neutrino mass. As an explicit example we have taken a symmetry with 
the charge assignments given in Table 1. The superpotential hence is then reduced to 

W = y3N5MH, + y[,N'5i,H-^ + y'^N'5^,H5 + 
5^/ 5^ + m^N'N 

yd 5m W H-^ + yu 10 10 H5 - fiH-^H5. (10) 

The superpotential in Eq. (5) is obtained identifying y3=Yj\f, y'i=Ki and y2=K2. Note that the term 
^H^H^ breaks the discrete symmetry explicitly (see discussion in Refs. [65, 71]). Furthermore, a somewhat 
problematic feature of this approach is, that the ^-problem of the MSSM is enlarged as not only /i, but 
also and rriM have to be few hundred GeV up to at most a few TeV. 

Usually, one argues that the /i problem can be evaded if one enlarges the MSSM Higgs sector by a 
singlet S. This leads to the NMSSM, where ^ = \{vs) and A is the coupling between the MSSM Higgs 
doublets and S. Indeed, one could in principle also obtain and m^v in this way in the desired order 
of magnitude. However, as we will discuss in the next section, this contradicts the requirement that the 
Weinberg operator should be forbidden. 



3 Implementation in the NMSSM, and the origin of a non-zero ^13? 

In the NMSSM, one introduces an extra gauge singlet S which couples to the Higgs doublets such that 
= \{S). Similarly one can generate in this way also the mass terms for additional particles in the model 
discussed above. For this, the generalization of the superpotential in Eq. (9) reads as 

W = yiN5^H^ + y2N5^,H5 + y3N5MH5 + 

y[ , N' 5^ + y'2 N' 5e H5 + y'3 N' 5m ^5 + 

A^/ S5m5^ + A^ 5% 5^ + XnSN'N + XnnSNN + Xn>n>SN'N' + 

yd 5m 10 % + y'd 5^' 10 % + Uu 10 10 H5 . (11) 

After electroweak symmetry breaking, the heavy mass scale is automatically set by {S) , which is of order 
of TeV. This has immediate consequences for the effective neutrino mass operators, which have now the 
schematic structure 

——LLHuHu, ——r{LLHuHd){HuHd). 

{S) {sy^ 

As a consequence {S) breaks any discrete symmetry under which it is charged. This prevents one from 
choosing a simple discrete symmetry group to avoid the d = 5 contribution to neutrino mass. Especially 
the term A^/ S 5m 5g , which leads to a mixing between the light and heavy mediators is problematic in 
this context, as it induces the Weinberg operator. 

We will prove now that one cannot forbid this term by a discrete symmetry that fulfills the required 
conditions. For this we start with three (yet) unconstrained charges as parameters 

qs = s, qHr, = 2h', QN = n. (12) 
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Multiplet 



Z3 charge 



5m H5 ifg N N' 5g 5g/ 10 S S' 



1 1 1 1 2 1 1 



Table 2: Charges for the fields of the model defined in Eq. (21). 



From the absolutely necessary terms in the superpotential we derive 

(SH^H-^) q^_=-s-2h' (13a) 

{SNN') q^, = -s-n (13b) 

{N5mH5) qM = -2h' - n (13c) 

(BmIO^s) gio = 4/1' + n + s . (13d) 

From the term (lOlOi^s) we obtain 

n=-s-5h' . (14) 

which leads to the following set of equations 

QHs = -s- 2h' , qn' = 5/i' , QM = 3h' + s , qio = -h' . (15) 

As a consequence we derive for the charges of the doublets 

(N'H^b^) q^ = -3h' + s (16) 

(7V'5g/F5) q^, = 7h' (17) 

(18) 



leading to 



qiSbMh) = 3s , (19) 



but we know that 3s = since we need a term in the NMSSM. This implies that one cannot forbid 
this unwanted term. Note, that this holds for every Abelian discrete symmetry group. As a consequence 
one can show along the same lines that every higher dimension operator for neutrino masses 



1 

{sy 



= ,^.^l^k+S LLHM{HuHdfS' (20) 



has to have the same charge as the Weinberg operator. We have also checked that similar problems 
appear if the Abelian symmetry is chosen as a product of two different cyclic groups ® Zjqi . 

The next possibility is the addition of a second singlet S' but this does not improve the situation either 
as one has to have interactions such that also S' obtains a VEV. However, from the term SSS' we obtain 
immediately g^/ = qs- As a consequence we find that the symmetry must be broken, which implies that 
also the Weinberg operator can be induced. However, if this breaking is soft, then the resulting Weinberg 
operator contributes only sub-dominantly to the neutrino masses and mixing angles. 

We consider the above model extended by a S' to exemplify the main points. The charges of the 
particles for the case of an unbroken symmetry are given in Table 2 and the superpotential^ at the SU{5) 
level is given by 



W = ysNbAi H5 + y[, N' 5^ % + y'2 N' 5^/ H5 + S' 5^/ 5^ + XnS' N'N 
+ yd 5m 10 % + yu 10 10 H5 + XsSH^H^ + kS^ + k'S'^ 

+ X'sS'H^H^ + y's N' 5m + y'd5^W + ■ ■ ■ , (21) 



^Note that this charge assignment would allow for a quadratic mass term for S' . However, we only allow for trilinear 
terms, which can be easily achieved, e.g., by adding an additional Z3 under which all fields have the same charge. 
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where the interactions in the first two lines respect the discrete symmetry and the ones of the last line 
break it. The X'gS' H^H^ term is a possibility to obtain a VEV for S' such that the mediators for the 
neutrino mass operators have their masses proportional to {S'). The second term of this line yields the 
Weinberg operator after integrating out ^ and 



=5 _ ysy's'^l 



{S') 



We obtain from the requirement 



'/\3 



(22) 



(23) 



that 7/3 < 10^® assuming (symmetry conserving) couplings of the order 10~^ and {S') ~ TeV. Hence we 
require all soft breaking couplings to be of order 10~® or smaller. We have checked that (5") is of O(TeV) 
in spite of X'^ being so small. 

It is interesting to observe that, for the soft coupling close to this upper limit, the neutrino mass 
matrix in Eq. (23) could receive a new (and potentially large) contribution from the d = 5 operator in 
Eq. (22). Suppose that the d = 7 operator has a tri-bimaximal flavor structure, where (for the sake 
of simplicity) the charged lepton mass matrix is assumed to be diagonal. Then the deviation from tri- 
bimaximal mixing can originate in the soft breaking terms in Eq. (21). We illustrate this effect with 
an exemplary flavor structure, not discussing the details of the underlying flavor symmetry. For reason 
of simplicity, we use a minimal scenario where the mass of the lightest neutrino is zero. For the d = 7 
contribution, we assume a flavor structure where we have three generations of the singlet fields and 
A'^' and the following structure of the coupling constants ^ 



I 

y\ 













, 2/2 = 2/2 




, 2/3 = 2/3! 











273 1/^/3 
l/y/Q 1/^/3 -l/y/2 
l/y/Q 1/^/3 l/y/2 



(24) 



where yi, y2 and are numerical (scalar) parameters, p = 7713/7712 and the mass matrices for N and 
A^' are diagonal. This choice of parameters will generate a tri-bimaximal mass matrix for the neutrinos 
at d = 7 



m 



d=7 



■ 2/3 



y'i(2/2) +2/2(2/'i 



(5)3 

Utb ■ diag(0, ?7i2, -ms) • U^^ 



yj 



(25) 



where Utb is the tri-bimaximal mixing matrix. We can write the total neutrino mixing matrix including 
small linear deviations from tri-bimaximal mixing by using the parametrization from Ref. [89] (note the 
different sign convention) 



C^PMNS 



i(l + s_a-r) i(l-is-a + r) -^(1 + a) 



(26) 



where sin0i3 = r/^/2, sin0i2 = l/\/3(l + s), and sin 6*23 = l/\/2(l + a). Since we want to generate the 
corrections to the neutrino mass matrix by the d = 5 contribution, we can write: 



mf-'^ + mfr^ , or 



(27) 



f^PMNS • diag(0, 7712, -7773) • C/pMNS = f^TB • diag(0,77i2, -7773) • C/tb + (2/3 (2/3)"^ + 2/3(2/3)""") (28) 



'Here we assume a similar fiavor structure as in Ref. [71] but three instead of two generations of A'^ and A'^'. 
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Solving this equation for the soft breaking couphngs 2/3, we obtain 



2/3 - % 



V6 , 

s — (a+r)p 
76 



-(2a+-r-s)-(2a-r)p^ 
2^3 

(2a-r-s) + (2a-r)p^ 
273 



(2a-r)-3r p^ \ 

373 * 

(2a-r)+3ap^ 

3^2 
(2a-r)+3ap^ 

373 



(29) 



to first order in r, s, and a. If we want to consider only corrections to ^13 {i.e., a = 0, s = 0), we have, 



( 

1 

'V6 




1 

"373 



/o 





3^2 ^ 

1 1 

2v^ 3\/2 







(30) 



with 0(^3) = 10 ^. This will give a correction to the tri-bimaximal mass matrix from the d = 7 operator 



that leads to the observed neutrino mixing with 
originate from a fiavor symmetry. 



?13 > 0. In a specific model, such a structure may 



4 Phenomenology of the additional d-quarks 

The promotion of the additional leptonic SU(2) doublets ^ and ^' to 5-plets 5^ and 5^/ of SU(5) implies the 
existence of additional d-quarks at the TeV scale as one wants to maintain gauge coupling unification. We 
will denote the corresponding mass eigenstates by D' , which is composed of d^' and d" , and L' which is 
composed of ^ and At the GUT-scale L' and D' have a common mass m^^GUT- Renormalization group 
efi^ects will lead to shift of the L' mass and the D' mass m£)i. At one- loop order the corresponding 
RGEs read as 



dt 



rrik 



-rrik 



(31) 



where the index i runs over the different gauge groups, t = ln(QVM2^2^), cf = (-4/15,0,-16/3) and 
c| = (-3/5,-3,0) for i = U{1)y, SU{2)l, SU{ci)c- Here we have assumed that one can neglect the 
Yukawa interactions compared to the gauge interactions. These equations can easily be solved and one 
obtains 



T^i,GUT 



(32) 



with hi = (38/5,2, —2). Clearly the precise ratio of these two masses at the electroweak scale depends 
on the particle content between the electroweak scale and the GUT scale as this changes the bi. In our 
model we obtain rriDi/m^ ~ 5 at Q = 1 TeV. Provided that these leptons can be observed at LHC up 
to mass of about 800 GeV [71, 90] the D' can have a mass of up to 4 TeV whereas LHC might observe 
them up to a mass of about 3 TeV see e.g. [90] and Refs. therein. 

Even though these quarks are heavier than the additional leptons, they will be stable in the model 
discussed so far as these particles are protected by the same symmetry that forbids the d = 5 operator. 
Stable heavy d-quarks, however, cause confiicts with cosmological constraints. Their presence during Big 
Bang Nucleosynthesis (BBN) would alter the observed abundances of the light elements in the universe 
(see, e.g., Ref. [91] for a review). Further bounds come, e.g., from direct heavy element searches in 
water [92]. Although these additional quarks can annihilate via gluons into SM quarks and gluons, it 
turns out this mechanism is not sufficient to lower the yield below the experimental requirements [93, 94] 
if their masses are below 2.5 TeV, e.g. the range interesting for LHC. To summarize, heavy stable d- 
quarks must have life-times much smaller than the age of the universe to avoid constraints from direct 
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searches, and annihilation processes that are efficient enough to lower their abundance below the bounds 
from BBN. These bounds do not affect particles that decay before BBN (r <C Is). 

In order to let the D' decay fast enough, one therefore needs additional soft breaking terms. For 
example the interactions of Eq. (21), discussed in the context of the NMSSM, lead to the two-body 
decays 

D'^H-u,H^d. (33) 

which result in life-times as small as 10^^*^ s to 10~^^ s, depending on m^)/, if the soft breaking couplings 
are of order 10~^. The signal for such a small life-time at the LHC is a displaced vertex. Measuring the 
corresponding decay length immediately gives the size of the corresponding coupling. 

One also has to take into account higher dimensional operators obtained from integrating out the 
additional gauge and Higgs fields at the GUT scale. For example, the decay of the heavy quarks via 
leptoquarks is possible and less suppressed than proton decay, due to the higher mass of the involved 
particles as the life time scales like rrv' . However, taking into account the current bounds on proton decay 
we find a lower bound on the life time of the heavy quarks that is still larger than the age of the universe. 
In addition there are the following two operators: 

e"e^^'""(5c)a(5M)fecde(^5).(%),Mm (34a) 
e'^^'^{b^UH^)j{H^)kmi^. (34b) 

Such operators can be obtained for example in the context of breaking Eq to SU{b) via an intermediate 
S'O(IO) with the help of an additional 5*0(10) 10-plet stemming from the 27-plet of Eq [95]. Neglecting 
the terms containing the heavy colored Higgs components and SUSY partners, operator (34a) leads to 
the terms 

(i7M«d, + e-L). (35) 

The second term gives a contribution to the Weinberg operator which however is only a subdominat 
contribution even if the couplings are of the size of a typical gauge coupling because it is suppressed at 
least by Mqut- The second operator (34b) leads to the term 

H^H+dflui. (36) 

Inserting the Higgs VEV, we obtain again the two-body decays of Eq. (33). However, even assuming that 
the couplings at the GUT scale have the strength of a gauge coupling, we find a lifetime of about 10^ s 
as the effective couplings are ■ Vd/McuT and ■ Vu/Mgut- Therefore, they give only small corrections 
to the soft breaking terms discussed above. 

5 Systematic review of the decompositions, and their SU(5) completions 

In the previous sections, we have discussed a specific realization of a supersymmetric model where the 
dominant contribution to neutrino masses is given by a d = 7 operator. The complete list of d = 7 
operators including the various mediators is given in [71], where the previous model is listed as #17. For 
convenience of the reader, we display here the possibilities in Table 3 where we also include the SU{5) 
multiplets containing the mediators. 

The corresponding SU{5) multiplets are either gauge singlets, 5-plets, 15-plets, 24-plets or 40-plets. 
The corresponding particles of the complete multiplet have to have masses at the TeV scale as not to 
destroy gauge coupling unification but only modifying the running of the gauge couplings. Using for 
example Ref. [96] one can easily calculate the additional contributions of these multiplets to the beta 
functions which we have summarized in Table 4. These values have to be compared with the MSSM 
beta functions bi = (33/5, 1, —3) for U{1)y, SU{2)l and S'C/(3)c, respectively. One sees that all but the 
5-plets give rather large contributions and one might wonder if such large contributions are consistent 
with the requirement of perturbativity up to the GUT scale. 
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# Operator Mediators SU(5) multiplets 

2 {H^iT^fL-){H^iT^L){HdiT^TH^) 3^, 3^ 1^, 1^, Sg 24, 24, (1), (1), 24 

3 {HjT^TL-)iHuiT^TL){Hdir^H^) 3^, 3^, Ig 24,24,1 

4 {-ie''''^)^uiT^T^L-){HjT^T''L){H,iiT^T^Hu) 3«, 3^, 3g 24,24,24 

5 {L-iT^TL){HdiT^Hu)iHuiT^THu} 3^1, 3^^, Ig 15,15,1 

6 {-i€abc)(l^iT^aL){HdiT^nHu){HuiT^TeH^) 3^i,3^i,3g 15,15^4 

7 iHJr^L-)iIAT^THd)iHjT^TH^) 1«, 1^, 3^i, 3^i, 3^1 1,1,15,1M5 

8 (-ie«^-)(if„ir2T'^L-)(LiTVF<j)(F„iT2r^F„) 3^ , 3^, 3«i, 3^^, 3^i 24,24,15,15,15 

9 {HjT^L-){iT^H^){L){HdiT^H^) 1^, 1^, 2^^/^, 2^i/2, Ig 1,1,5,5^1 

10 {HJt^tL'=) {it^tHu) (L) {HdiT^Hu) 3^, 3^, 2^,^,, 2^,^,, Ig 24,24,5,5,1 

11 {HuiT^T7-){iT^Hu){fL){HdiT^THu) 1^, 1^, 2«,/2: 2^1/2. 3g 1,1,5,5^24 

12 {HjT^T'^L-){iT^T''Hu){T''L){HaiT^T'>Hu) 3«, 3^, 2«^/2, 2f:,/2, 3g 24,24,5,5,24 

13 {HuiT^L-){L)(iT^Hu){HdiT^Hu) 1^, 1,^ 2;,/2, Ig 1,1,5,1 

14 (HuiT^TL-){fL){iT^Hu){HdiT^Hu) 3^, 3^, 2;i/2, Ig 24,24,5,1 

15 {HuiT^L-){L) (ir (grfir l^s 1^, 2;,^,, 3g 1,1,5,24 

16 {HjT^T''T?){T''L){iT^T'>H^)iHdiT^T^Hu) 3^, 3,f , 2;^/^, 3g 24,24,5,_24 

17 {HuiT^L-){Hd){iT^Hu){HuiT^L) 1^^.1,^.2^^/2,2^/2 1,1,5,5 

18 (HjT^TL-){fHd){iT^H,,,)(HjT^L) 3^, 3^, 2«i/2, 2^i/2, 1^', 1^ 24, 24, 5, 5, (1), (1) 

19 {HjT^L-){Hd){iT^fH^){H^iT^fL) 1^,1^,2^^/2,2^^/2,3^,3^ (1), (1), 5, 5, 2_4, 24 

20 {HjT^T''L-)iT''Hd){iT^T''H^){HjT^T^L) 3g^, 3^, 2g^/2, 2^:^/2, 24,24,5,5 

21 (L^iT^T''L){HjT^T-){T'Hd){HuiT^T'Hu) 3^1, 2;,/2 , 3^1 15,5,15 

22 {L-iT^"L){HaiT^T''){T^Hu){HjT^T^H^) 3^1,2^3/2,3^1 15,40,15 

23 _{L-iT^TL){HuiT^f){Hu){HdiT^R,) 3^+i, 2;^/2, Ig 15,5,1 

24 (L-irV''L)(i?„irV)(T''i7„)(ffdirVi7„) 3^i, 2;^/2, 3g 15, 5_, 24 

25 {HaiT^Hu){L-iT^)ifL){HuiT^THu) Ig, 2^^/2, 2g^/2, 3^^ 1,5,5,15 

26 {HdiT^T-H^(l?iT^T-){T'L){HjT^T'Hu) 3g, 2^,/2, 2^,/2, 3^1 24,5,5_,15 

27 (if„iT2L-)(ir2Fd)(fL)(if„iT2fif„) 1^, 1^, 2^,/2, 2^^/2, 3^1 1,1,5,5^15 

28 iHjT^T-L^{iT'T-Hd){TbL){HjT'T'H^) 3«, 3^, 2«^/2, 2^^/2, 3^1 24, 24, 5, 5, 15 

29 {HuiT^L-){L){iT^THd){HuiT^TH^) 1^,1^,2^1/2,3^1 1,1,5,15 

30 (g„irV°L^) (r°L) (ir V'ffd) (gjr V'g^) 3g^, 3^, 2;,/2, 3^i 24,24,5,15 

31 {mT^T''Hd){iT^T''Hu)iT''L)iH^iT^T''Hu) 3^1,3^1,2^^/2,2^^/2,3^1 15,15,5,5,15 

32 (i-irV«ifrf)(T«i)(iTV''ff„)(if„iTV*'ifJ 3^:i, 3fi, 2;3/2, 3^1 15j5,40,15 

33 _(L=iT2fJ?d)(iTVi?„)(F„)(7J„iT2L) 3^i, 3«i, 2^1/2, 2«i/2, 1^, 1^ 15^5, 5, 5, 1, 1 

34 (MrMifd)(irMiJ„)(r''£r„)(J?„iTVL) 3^i, 3^1, 2^:1/2, 2^1/2, 3^, 3^ 15,15,5,5,24,24 



Table 3: Decompositions of the d = 7 operator LLHuHuHaHy, at tree level. We use the following notation 
for the mediators: Xy. The X describes the SU(2) nature, i.e., singlet (1), doublet (2), or triplet (3). 
The superscript L denotes a left- (L) or right- (R) handed fermion, or a scalar (s). The subscript Y 
represents the hypercharge Y = Q — I^. The SU(5) singlets in parentheses can be already contained 
in the 24-plets. 
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Multiplet 


5 


15 24 40 




1/2 


7/2 5 11 



Table 4: Contributions of the various SU{5) multiplets to the MSSM beta functions. 



At the one-loop level the RGE of the gauge coupling have the well-known solution 

«G = T72- , (37) 

where the Q is the scale at which the additional particles are included, Mq ~ 2 • 10^^ GeV is the GUT 
scale and ac the value of the gauge couplings at Mq- Taking only the MSSM one gets ac — 1/25. 
Adding a contribution Abi to bi in Eq. (37) we find Abi = 5 already implies that ac = 1, see also [97] 
for a related discussion. As the 5-plet and the 15-plet are complex representations, they have to come 
always paired with a 5 and 15, respectively, as otherwise one cannot put a gauge invariant mass term in 
the superpotential. Inspecting Table 4 we thus conclude that we can add at most five 5, 5 pairs or one 
24-plet if one is willing to accept oq = 1. If one would require ac < 1/2 then one is restricted to at most 
four 5, 5 pairs. In practice the number is smaller as two- loop effects increase the value of ac [84]. As a 
consequence, from the possible operators in Table 3 only operators #1, #9, #13 and #17 are consistent 
with perturbativity up to the GUT scale. We note for completeness that in case of operator #9 one has 
to include a corresponding 5 as otherwise no mass term can be added to the MSSM superpotential. Also 
for operators #9 and #13 the discussion as for operator #17 in the previous sections can be repeated 
resulting in essentially the same findings. Note that decompositions #1, #9, and #13 contain a singlet 
scalar in the decomposition. Looking at the according superpotentials 

W#i = 5m10H-^ + 10WH5 + SSS +S H^H-^ + H^BmN + H^H-^cP + 

N'N'cp + N N'S +^(pS (38a) 

W#3 = 5mWH-^ + WWH5 + SSS +S + //55Af24, + H^H-^cj) + 

24b2ib(p + 2ia2ibS +<I)(I)S (38b) 
W#i3 = 5mW H-^ + WW H5 + S S S + S H^H-^ + SmH^N + N'5m5s+ 

5sH5(f) + H^H-^cf) +NN'S + 5s5sS +(I)^S , (38c) 

where the (j) are the scalar singlets and for the other fields we use the notation introduced in the previous 
sections. Note, that we only give the superfields but not the details of the couplings. All three super- 
potentials contain the term H^H^cj), which eventually induces a vev for the scalar component of (j). In 
these cases, a lower than d = 7 dimensional operator will dominate neutrino mass. This means that our 
example in Fig. 1 is the only tree level decomposition which is consistent with perturbativity up to the 
GUT scale and neutrino mass generation by a leading d = 7 contribution. 



6 Summary and conclusions 

We have discussed neutrino mass generation from a d = 7 effective operator, while the d = 5 Weinberg 
operator is suppressed due to a discrete symmetry. The additional suppression by u^/A^p compared to 
the Weinberg operator reduces the new physics scale Anp required to describe the smallness of neutrino 
mass, potentially down to the TeV scale, which makes the mediators accessible at the LHC. We have 
especially focused on TeV completions in the MSSM and NMSSM frameworks. In order to achieve gauge 
coupling unification, we have embedded such models in an SU(5) GUT. 

We have chosen one specific example for the decomposition of the d = 7 operator in Eq. (4) , see Fig. 1 , 
which has been discussed earlier in the context of its LHC phenomenology [71]. It involves two fermionic 
SU(2) singlets, and two fermionic SU(2) doublets. In this case, not only the mediators may be probed at 
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the LHC, but also sensitivity to lepton flavor and perhaps even lepton number violation can be achieved. 
In this study, we have demonstrated that this example is peculiar for two reasons: 1) it can be completed 
in an SU(5) and described perturbatively up to the GUT scale, and 2) it does not contain any singlet 
scalars which take a vcv and thereby generate neutrino mass by a lower dimensional operator. It is, in 
fact, the only example which satisfies these two requirements. 

The embedding of the fermionic doublets into SU(2) multiplets implies that additional heavy d-quarks 
are to be introduced, which are strongly constrained by early universe cosmology, such as big bang 
nucleosynthesis. Therefore, fast enough decay mechanisms for these d-quarks need to be identified. The 
decay via leptoquarks or higher dimensional operators, which may be generated from physics beyond 
the GUT scale, is possible, but nevertheless too slow to satisfy the current constraints. The possibility 
in agreement with cosmological bounds is a decay via a soft (discrete) symmetry breaking term, which 
could be large enough to allow for a sufficient decay rate. 

For the specific case of the NMSSM, we have demonstrated that the NMSSM scalar vev {S) breaks 
the discrete symmetry spontaneously, which will lead to a (dominating) neutrino mass term from the 
Weinberg operator. A simple way out has been shown to introduce an additional scalar and soft breaking 
terms, which can be small enough such that the d = 5 contribution to neutrino mass is suppressed. 
We have also shown that a soft breaking term may describe the origin of deviations from tri-bimaximal 
mixing, and this may be the origin of the non- vanishing magnitude of ^13. 

In conclusion, we have discussed a model which described neutrino mass from physics at the TeV 
scale, and which is potentially LHC-testable. We have embedded the model in an SU(5) SUSY-GUT, 
which means that wc have confirmed gauge coupling unification and pcrturbativity up to the GUT scale. 
Finally, we have established that this model requires the fast enough decay of heavy d-quarks, which are 
needed for the GUT completion of the theory. Therefore, interesting constraints will not only be obtained 
from LHC, but also from early universe cosmology. 
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